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Abstract 

We present a model of polymer growth and diffusion with frustration mechanisms for density 
increase and with diffusion rates of Arrhenius form with mass-dependent energy barriers T{m) ~ 
(m — l) 7 . It shows non-universal logarithmic coarsening involving the exponent 7. Strong-glass 
behavior is found in the typical times for disappearance of all polymers up to a given length, 
without reference to the equilibrium states of the macroscopic system. These features are predicted 
by numerical simulations, scaling theories and an analytic solution of the master equation within 
an independent interval approximation, which also provides the cluster size distribution. 
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I. INTRODUCTION 



One of the most remarkable features of glassy systems is the rapid increase of the re- 
laxation time r to equilibrium states as the temperature T is lowered. In the most simple 
glasses, called strong glasses, an Ahrrenius behavior r ~ exp {A/T) is observed, while in 
fragile glasses more complex temperature-dependences are found. This slow relaxation is 
expected to be accompanied by a slow growth of correlated domains. Several microscopic 
models have already been proposed to represent such features jll Q] • For instance, strong 



glass 
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jehavior was found in the spin-facilitated models introduced by Fredrickson and Ander 



4 and several glassy features were found in models with kinetic constraints 
in which initial random systems evolve to equilibrium states via the slow diffusion of deffects. 
On the other hand, fragile glass relaxation to equilibrium states with r ~ exp (B /T 2 )) was 
recently shown in a spin chain with asymmetric kinetic constraints 0, 1(| . The dynamic rules 
assumed during the non-equilibrium evolution of these systems are derived from statistical 
equilibrium conditions. For instance, in the spin-facilitated models, the final concentration 
of defects is given by the canonical distribution exp (-E/T) at low temperatures and this 
is the origin of the Arrhenius form of the relaxation time in those systems. 

An alternative scenario for the onset of anomalous coarsening and glassy behavior is 
suggested in this paper with the analysis of a polymer growth model in one dimension. 
The slow dynamics in this model is a consequence of the interplay between slow activated 
(Arrhenius) diffusion of clusters and frustration of density increase. Cluster diffusion occurs 
in thermal contact with the surroundings, with energy barriers increasing with cluster length. 
Density increase is represented by the deposition of new particles in the line, but it is not 
allowed at small vacancies between the clusters. Thus, the dynamic rules of this model 
do not make any reference to equilibrium macroscopic states and the physical motivation 
of those processes contrasts with the somewhat artificial stochastic rules of other simple 
models with similar dramatic slowing. 

The rules of the model, illustrated in Figs, la and lb, prescribe as follows the influ- 
ence of polymer length on activated diffusion, the suppression of density increase, and the 
irreversible polymer aggregation. A cluster (polymer) can move one lattice spacing to the 
right or to the left with diffusion rate given by r = exp(—T/T), where the energy barrier 
r increases with polymer length and T is the temperature (Fig. la). We will assume that 
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r ~ (m — l) 7 , where m is the polymer length (mass) and 7 > 0. Deposition of hard core 
particles, which represents density increase, is allowed only at sites with one empty nearest 
neighbor, with deposition rate F — 1 (Fig. lb). Aggregation of a particle to a cluster and 
of two clusters is irreversible and occurs upon any contact of nearest neighbors. Thus, after 
a diffusion event, the most typical situation is the formation of a larger cluster (Fig. la). 
This represents an ideal polymerization process, with no energy barrier for the formation 
of a new bond between neighboring particles and an infinite barrier for the reverse process. 
The linear increase of T with m (7 = 1) would correspond to polymers stretched along a 
surface, with V being the sum of adsorption energies of all monomers. Although not directly 
related to the present problem, it is also relevant to recall that energy barriers with 7 m 0.5 
are found in desorption of linear alkanes from graphite surfaces jll|. Consequently, our 
numerical study will be concentrated in systems with 7 ~ 1. 

(&) r=exp[-(m-l) Y /T] 
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FIG. 1: (a) Diffusion of a polymer of length (mass) m, with the corresponding rate r, and the two 
possible final configurations after the polymer move, (b) Allowed deposition processes at vacancies 
with a neighboring vacant site, with the corresponding rate, and the forbidden deposition process, 
in which the vacancy has two occupied neighbors. 

The restriction of the model to a low- dimensional structure helps to obtain the analytical 
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and numerical solutions. This is certainly a simplification for many systems, in which the 
dimensionality may play an important role. However, glassy behavior has been seen in 
several low- dimensional systems, such as polymer films |l2j , and the formation of clusters of 
fast-moving particles with dimensions near or below 2 was observed in colloidal glasses [13| . 
Consequently, we believe that the qualitative scenario introduced here may find applications 
to real systems. 

Using scaling arguments and numerical methods, we will show that the model presents 
non-universal logarithmic coarsening, in which the average cluster (polymer) length slowly 
increases as (m) ~ (Tlnt) 1 ^ 7 , starting from an initial random low-density configuration. The 
model also presents strong-glass features of the characteristic times for eliminating structures 
of a given length: after a time r ~ exp [<P /T], all clusters of sizes of order d or smaller will 
disappear. Consequently, the anomalous coarsening and the Arrhenius-type relaxation are 
observed during the non-equilibrium system evolution, which involves frustration of density 
increase and thermally activated diffusion at the microscopic level. We will also solve the 
master equation of the process within an independent interval approximation, using the same 
map ping onto a column picture adopted in the study of related systems in Refs. and 
|lq . However, the treatment of cluster size distributions of the present model has novel and 
less trivial aspects. It will confirm the above results and provide cluster length distributions 
that remarkably differ from those related systems and that, as far as we know, were not 
previously obtained in other systems with logarithmic coarsening. 

The rest of this paper is organized as follows. In Sec. II we use scaling arguments to 
predict the logarithmic coarsening and the glassy behavior and confirm these results with 
numerical simulations. In Sec. Ill we solve the master equation within an independent 
interval approximation, providing the cluster length distributions. In Sec. IV we summarize 
our results and conclusions. 

II. SCALING OF CLUSTER LENGTH AND GLASSY BEHAVIOR 



Scaling arguments can be used to predict the slow coarsening of this mode 



same lines of Ref. 
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Q, which were _ly applied to magnetic systems QQandto 
related non-equilibrium models [I4J,115||- For simplicity, we refer to the average cluster mass 
as m. 
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FIG. 2: Illustration of the coarsening dynamics for the construction of the scaling theory: (a) 
three typical neighboring clusters are shown, with lengths of the order of the average cluster size 
m; (b) cluster B moves to the left, forming a larger cluster after aggregation to cluster A, and a 
new particle is deposited at one of the sites of the double vacancy between A + B and C; (c) the 
final configuration after the deposition of a new particle at the left site of the double vacancy. 

In Fig. 2a, we show a configuration with clusters of lengths typically of order m, named 
A, B and C, with single empty sites between them. The time necessary for cluster B to 
move is of order At = exp (m 7 /T). If it moves to the left (Fig. 2b), then clusters A and B 
coalesce. Since the diffusion rates of clusters A + B and C are very small for large m, while 
the deposition rate is F — 1, a new particle will immediately be deposited in one of the sites 
of the double vacancy (Fig. 2b). One possible configuration after the deposition of a new 
particle is shown in Fig. 2c. From the initial to the final configuration (Fig. 2a to Fig. 2c), 
the average cluster length increase was of order m. Thus we obtain 

dm Am m 

~dt At ~ exp (m-y/T)' ^ ' 

We will be mainly interested in the long-time regime in which cluster diffusion is suffi- 
ciently slow, even for small 7. In this case, the exponential factor at the right-hand side 
of Eq. is much larger than m, which may be neglected at first approximation. Thus, 
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integrating Eq. (JTJ, we are led to the scaling of the average cluster length as 

m~(T In tf h . (2) 

The leading correction due to the neglected term in Eq. ((TJ is proportional to 
In (lnt)/(lnt) 1-1 ^ 7 . However, it is not expected to be the true leading correction to the 
dominant scaling of the model because Eq. (P) is itself an approximation, which omitted 
further corrections. 

Numerical simulations of this model confirm the logarithmic coarsening. In Fig. 3 we 
show the time evolution of the average cluster length for three sets of values of 7 and T, in 
which that scaling is observed in up to 7 decades of time. The data for 7 = 1, with T = 1 
and T = 2, are averages over 3 x 10 4 different realizations in lattices of length L = 10 4 . The 
data for 7 = 0.5, in which much larger clusters appear at small times, are averages over 200 
realizations in lattices with L = 5 x 10 4 , up to time t = 10 5 . 
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FIG. 3: (m) 7 versus log (t) for: 7 = 1 and T = 1 (squares); 7 = 1 and T = 2 (triangles); 7 = 0.5 
and T = 1 (crosses). 

In Fig. 4 we show the data for the same values of T and 7 scaled in the form (m) / (T In t) 1 ^ 7 
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versus 1/ Int. As t — > oo (1/ hit — * 0), Fig. 4 suggests an universal amplitude in the scaling 
relation (j2J), although corrections to the dominant scaling clearly depend on 7 and T. 
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FIG. 4: Scaled average cluster size as a function of inverse of lni for: 7 = 1 and T = 1 (squares); 
7 = 1 and T = 2 (triangles); 7 = 0.5 and T = 1 (crosses). 

The logarithmic coarsening leads to a time evolution of the density with a logarithmic 
correction, p = 1 — 1 / (T In t) . This type of time-dependence was also found in experiments 



of compaction of granular systems [19] and theoretical models to describe them |20|,|2l|, but 
temperature played no role there. 

Another remarkable feature of this system seen in simulations is the absence of clusters 
whose reduced lengths 

y = m/{T In t) 1/7 (3) 

are smaller than the most probable value y = a (a —>■ const for large t). In order to 
understand this feature, consider a set of clusters with reduced size y = a — e, with any 
e > 0, whose diffusion rate is t^ a ~^ 1 . This value is larger than the diffusion rate of the 
clusters with y = a by a factor t n ^ a . Thus, at large t, compared to the most probable 
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clusters, the clusters with y = a — e are much more mobile and are present in smaller 
number in the line. So they will merge into the slower larger neighboring clusters almost 
instantaneously within the typical time scale for diffusion of the most probable clusters 
(length y = a), which causes the removal of clusters with y = a — e. 

This is consistent with the simulation results shown in Fig. 5, in which we plotted the 
reduced probability of cluster length m, (Tlnt) 1 ^ 7 P (m), as a function of the reduced length 
y, for 7 = 1 (two temperatures in the main plot) and 7 = 0.5 (in the inset). Below a certain 
value y = a « 1, the corresponding probabilities rapidly decrease to zero and the cluster 
size distribution is asymptotically discontinuous at that point. A discontinuity in the slope 
of the distribution at y ~ 2a is also suggested in Fig. 5 and will be rigorously justified in 
Sec. III. Further slope discontinuities will be predicted, but they are not so clear in Fig. 5 
due to the larger fluctuations in the data for large y. 
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FIG. 5: Scaled probability of clusters of length (mass) m as a function of the scaled length y 
for 7 = 1, with T = 1 (solid curve) and T = 2 (squares), at t = 10 8 . The inset shows the same 
quantities for 7 = 0.5 and T = 1, at t = 10 (not shown in the main plot to avoid superposition of 
many data points). 
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From the above results, the characteristic time for clusters of a given mass m to be 
completely eliminated from the system is t — exp [(m/a) 7 /T]. The connection of the glassy 
behaviour to the thermally activated diffusion of microscopic structures is clear, in contrast 
to other kinetically constrained models whose dynamic rules are determined by the statistical 
equilibrium of the system. 



III. ANALYTICAL SOLUTION 



Now we turn to the analysis of the model starting from a version of the master equation, 
following the same lines as in the solution of related models in Refs. Q and Q- 

In our original problem, each site on the line is occupied by one or zero particles, but we 
are interested in the evolution of cluster length. Thus, the analysis of a master equation is 
more easily set up by reformulating the process using a column picture, in which a column 
of height m represents a cluster of size m together with its adjacent vacancy on the right. 
This mapping is illustrated in Fig. 6. 



(b) 




FIG. 6: Example of particle-hole configuration on a line (a) and the map (dashed arrows) into a 
column problem (b). Each cluster and the vacancy at its right side in (a) correspond to a column 
in (b) with the same mass. 



Fig. 6 clearly shows that the lattice length Lq in the column picture is smaller than the 
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length L in the original problem (particle picture): Lq = L — M, where M is the total mass, 
for periodic boundaries. As time increases and new particles are deposited, L decreases. 
However, new particles are deposited only at double vacancies of the particle picture (Fig. 
lb) and the latter correspond to single vacancies in the column picture (Fig. 6). Thus, the 
decrease of L (or mass increase) is related to the probability Pt(0) of an empty site in the 
column problem as Lo{t+ ^ t ) 0{t) = ~P(0) [2 - P(0)]. 

The probability Pt{m) that a randomly chosen cluster in the column picture has size m at 
time t is given in terms of the clusters numbers, N(m,t), by Pt{m) = N ^t) • g a i n /l° ss 
from in/out processes provides a master equation for clusters numbers as 

N(m,t + 1)-N(m,t) = L {B m + E m ), (4) 

where B m comes from the deposition processes and E m comes from diffusion of clusters or 
isolated particles. It can be rewritten in terms of cluster length probabilities as 

(1 - P(0)) 2 P m (m) - P t (m) = B m + E m , (5) 

which allows for the changes in L . 

In an independent interval approximation (IIA) in which joint probabilities are factorized, 
the contribution from deposition is given by 

B m = P(0) [29 (m - 2) P t (m - 1) - 29 (m - 1) P t (m) + 5 m>1 P t (0) - 2S mfi ] . (6) 

Here, the factor Pt{0) accounts for the fact that deposition is allowed only at single vacancies 
in the column picture (Fig. 6). The first term inside brackets gives the rate of production 
of clusters of size m from deposition at one of the sides of a cluster of size m—1, for m > 2, 
while the third term gives the rate of production of a cluster with m = 1 in a double vacancy 
(triple vacancy in the particle picture, with the new particle being deposited at the central 
site). The second term corresponds to the loss of a cluster of mass m due to deposition at 
one of its sides and the last term the loss of vacancies due to deposition at one of the two 
corresponding sites in the particle picture. 

The contribution from diffusion in Eq. (jHJ) is given by 

m oo 

E m = £/(r)P t (r)P t (m-r) + <5 m , £/(r)P 4 (r) 

r=l r=l 

oo 

-6(m - l)f(m)P t (m) - P t (m) £ f(r)P t (r). (7) 

r=l 
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In Eq. (JJJ), the function f(m) = exp [— (m — 1) 7 /T] gives the diffusion rate for a cluster 
of mass m. The first term in Eq. (j7J) gives the rate of production of a cluster of mass m 
due to the coalescence of smaller clusters after diffusion of one of them. The second term 
corresponds to the creation of a single vacancy (double vacancy in the particle picture) due 
to cluster diffusion - see Fig. la. The loss terms correspond to coalescence of a cluster of 
mass m and a neighboring cluster: the third one accounts for diffusion of the cluster of mass 
m and the fourth one accounts for diffusion of the neighboring clusters. 

In the scaling limit of large t and large m (coarsening), the probability Pt{m) must have 
the form 

1 in 

p ^ m) = T{t) 9{y) ' y = WY (8) 

with some function L(t) characterizing the growing typical cluster size. 

Figs. 2a-c show that the survival time of a double vacancy is of order 1/F = 1, while a 
typical cluster move occurs on a much larger timescale, of order exp (m 7 /T). Consequently, 
the onset of a double vacancy (in the particle picture) is a rare event. Nevertheless, since the 
double vacancy mediates deposition, its probability is required in the equation for the one- 
variable scaling function g, and the same scaling arguments of Sec. II lead to the expected 
scaling of the double vacancy probability in the particle picture as P t (0) = (C/2) jt (see 
e.g. Ref. H). 

Thus we obtain an equation for the scaling function as 



T 



, d 9 
9 + V-r 
dy 



C<l Cdi!j) -h(y)-g(y)I(oo) + 5(y)I(oc) + 
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dy'g(y-y')h(y'), (9) 



where h(y) = g(y) exp [— (Ly) 1 /T] and I(oo) = J °° dy'g (y') exp [— (Ly 1 ) 1 /T}. The first term 
on the left-hand side (LHS) of Eq. (JHJ) follows from those in Eq. (jSJ in the long time, 
continuous limit. The second term on the LHS and the first one on the right-hand side 
(RHS) correspond to deposition (Eq. |HJ) and the other terms on the RHS correspond to 
cluster diffusion (Eq. |7j). 

Eq. and the correpon ding equation for the generating function contrast to those of 
the models of Refs. Q| and jlfij] because a simple power-counting in Eq. (jSJ) is not enough 
to provide the scaling of the average cluster size here. From the structure of Eq. and 
the forms of h(y) and I(oo), we expect that the balance of dominant terms will be possible 
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if 

j = exp[-[L(t)a]yT]a(L(t)), (10) 

where a is some constant (see below) and the function a accounts for possible power-law 
subdominant factors. Notice also that the dominant factors of Eq. fllOJ) are the same as 
those of Eq. (fT|) . obtained from simple scaling arguments. 

As t — > oo, we expect L(t) — > oo, but the last term in the RHS of Eq. will provide 
a divergent contribution for y' < a. This contribution is not balanced out by other terms 
unless g(y) = for y < a. With this assumption, we obtain a dominant term on the LHS 
of Eq. (JHJ) with a factor e~^ La ^ '/ T and, in the integral at the RHS, a factor e~^ Lyn>1 ^ T , which 
provides a dominant cancelling term for y > a. 

The solution of Eq. (fTT)|) is exactly the logarithmic growth of L(t) obtained from scaling 
arguments (Eq. |2J), thus providing the expected scaling of average cluster length. The above 
assumption for the function g confirms, within the IIA, the elimination of all clusters with 
masses smaller than aL{t) ~ (Tint) 1 / 7 . 

The scaling function has the form 

g(y)=e(y-a) gi (y), (11) 

which leads to an equation for the function g\{y) as 

gi+y^+Cgi-C6( y ) = ^2 [ 9l 9( y -a)- 9l (y - a)9(y - 2a) + 5{y - a) - 8{y)\ , (12) 

where k is constant. Consequently, a discontinuity in the slope of g(y) is obtained at y = 2a. 
Subsequently, further discontinuities appear at y — 3a, y = 4a and so on. The function g 
can be shown to be continuous except at y — a, where the discontinuity implies k = 1 in 
Eq. (|T2*|l . These findings agree with those in Fig. 5, in which the probabilities of cluster 
lengths with y < a ~ 1 tend to zero asymptotically and the slope of the scaled probability 
clearly changes at y ~ 2a. 

From Eq. (|T2*|) we obtain the distribution for a < x < 2a as 

g(y) = Ay^-^\ (13) 

where A is a constant. In Fig. 7 we show log P(m) versus logy for 7 = 1 and T = 1, at time 
t = 10 8 , with a linear fit that confirms the power-law decay predicted in Eq. (fT3*j) . Notice 
that, although the scaling arguments of Sec. II were capable of predicting several features of 
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the model, the analytical solution of the master equation, even within the IIA, is essential 
to predict the shape of the scaled cluster size distribution. 

-3 i 1 1 
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ln(y) 

FIG. 7: In [P (m)] versus In (y) for 7 = 1 and T = 1, at t = 10 s , and a linear fit of the data in the 
region 1.4 < y < 2.2. 

For 2a < y < 3a, the distribution is g(y) = By^ a ^ c+1 ^ _ A ( y _ a ^(fl)-(c+i)] ^ but the 
low accuracy of our data in that region does not allow a reliable test of this form. 

IV. CONCLUSION 

We presented a one-dimensional model of polymer growth and diffusion with frustration 
mechanisms for density increase and with diffusion rates of Arrhenius form, with mass- 
dependent energy barriers. The non-universal logarithmic coarsening involves the exponent 
7 and can be predicted by scaling arguments. The model also shows strong glass behavior 
of the characteristic times for elimination of all polymers up to a given length. These 
features are confirmed by the solution of the master equation within an independent interval 
approximation (IIA). This solution also provides the distribution of cluster sizes, which shows 
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discontinuities at a sequence of points and a power-law decay for sizes near the most probable 
one. These findings are confirmed by numerical simulations with good accuracy. 

The logarithmic coarsening and Arrhenius-type behavior were also obtained in other sys- 
tems, but our model presented a different scenario for the onset of these properties, in which 
there is no reference to equilibrium states. Instead, all these features are obtained during 
the non-equilibrium system evolution assuming a thermal contact between the clusters and 
their surroundings at a temperature T, but not an approach to an equilibrium macroscopic 
state with that temperature. 

The elimination of clusters with lengths below a certain threshold is also present in the 
paste-all model of Derrida et al (2^], although that is a direct consequence of the model 
prescription, namely the elimination of the smallest cluster of the system at any time. 
Instead, in the present model that elimination is obtained only in the scaling limit (long 
times, large clusters), in which diffusion coefficients of small clusters become infinitely larger 
than those of large clusters. The same type of distribution appears in the East model 0, [l0| , 
in which the rules for spin flips are asymmetric (an artificial but essential ingredient for 
obtaining slow dynamics there). On the other hand, a power-law coarsening and fragile 
glass behavior are obtained in the East model, which contrasts to our findings. 



The successful application of the IIA 



with particle detachment from clusters 
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or the present problem and for previous models 



15| also deserves some comments. In these mod- 



els, the coarsening arises from coalescence of clusters and the reverse process is increasingly 
improbable (altogether improbable in the present case). Since the initial state is without 
cluster-cluster correlations, no correlations between the masses/ lengths of neighbouring clus- 
ters can build up, so the IIA becomes exact in the late coarsening limit. This interpretation 
is corroborated by the successful comparison of IIA predictions and simulation results. It 
seems that IIA fails only when we focus on rare processes which occur in narrow time win- 
dows in which the system is dominated by reversible and highly correlated processes jlfil ]. 
which is not the case in the present model. 
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